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Abstract 

Domain Wall Fermions utilize an extra space time dimension to provide a method 
for restoring the regularization induced chiral symmetry breaking in lattice vector 
gauge theories even at finite lattice spacing. The breaking is restored at an expo- 
nential rate as the size of the extra dimension increases. Before this method can be 
used in dynamical simulations of lattice QCD, the dependence of the restoration rate 
to the other parameters of the theory and, in particular, the lattice spacing must be 
investigated. In this paper such an investigation is carried out in the context of the 
two flavor lattice Schwinger model. 
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1 Introduction 



When fermions are discretized on a d-dimensional lattice they "double" producing 2 d species 
for each flavor. In order to remove the unwanted degrees of freedom special care must be 
taken. For a vector theory, like QCD, two methods have been used to deal with this 
problem, but both break the global symmetries of the continuum theory. Wilson fermions 
[.fj are implemented by adding an irrelevant operator to the action. This operator makes all 
but one of the species heavy (with masses close to the cutoff). For the Nf flavor QCD this 
operator breaks the SU(Nj)l x SU(Nf)n chiral symmetry down to SU(Nf). This explicit 
breaking is severe and requires fine tuning of the bare quark mass in order to obtain a 
massless theory. Even then the size of the breaking is proportional to the lattice spacing and 
only close to the continuum limit the explicit breaking becomes small. Staggered fermions 
break the SU(Nf) L x SU(Nf) R chiral symmetry down to £7(1) x £7(1). Because of the 
remnant of chiral symmetry the massless theory can be reached by simply taking the bare 
quark mass to zero. However, the flavor symmetry of the theory has been compromised 
and is also only recovered as the continuum limit is approached. 

Despite these problems both methods have been very successful in describing the light 
hadron spectrum at zero temperature. However, both methods have difficulties in studying 
the finite temperature phase transition. Wilson fermions have a complicated phase diagram 
that, at the presently accessible lattice spacings, makes it hard to extract the relevant 
physics. Staggered fermions, because of the exact remnant of chiral symmetry, do not suffer 
from this problem. However, at the presently accessible lattice spacings, the breaking of 
flavor symmetry makes two of the three pions heavy. This can have important physical 
consequences since the transition temperature is of the order of the pion mass. For a review 
on the finite temperature phase transition with both types of fermions the reader is referred 
to [[| and references therein. 

A few years ago a new method for discretizing fermions was developed in order to address 
the more difficult problems associated with chiral gauge theories Q . In the following years 
this method was further developed (see || and references therein) with important progress 
in the development of chiral gauge theories 0,0. The basic idea follows from the fact that 
a massive vector theory in 2n + 1 dimensions, with a mass term that changes sign along 
the 2n + 1 dimension, develops a massless chiral zero mode that is exponentially bound 
along the 2n + 1 direction to the region where the mass changes sign. From the point of 
view of the 2n dimensional world this is a chiral fermion. This region is called Domain 
Wall and this type of fermion is called Domain Wall Fermion (DWF). When such a theory 
is discretized species doubling also occurs. However, since the 2n + 1 dimensional theory is 
vector-like the extra species can be removed with the addition of a standard Wilson term. 
The resulting theory has a single chiral fermion exponentially bound to the wall. If, for 
practical reasons, the 2n + 1 dimension is made finite with periodic boundary conditions 
for the mass then the mass must change sign one more time. In that region (anti-wall) 
an exponentially bound chiral zero mode with opposite chirality appears. As a result, the 
theory becomes vector-like. Different types of boundary conditions yield similar problems. 
In order to preserve the single chiral mode, the 2n + 1 dimension must be kept infinite. 
At first sight this may seem impractical. However, Narayanan and Neuberger developed a 
method, called the Overlap formalism, that makes it possible to deal with this infinity ||]. 
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The Overlap formalism develops a transfer matrix along the 2n + 1 dimension and an 
associated Hamiltonian. The gauge fields are defined only on the 2n dimensional space and 
are taken to be independent of the 2n + 1 coordinate M. In essence, the extra dimension is 
treated as a complicated flavor space. The resulting formalism involves two Hamiltonians, 
one for the region of positive mass and one for the region of negative mass. The chiral 
determinant is the determinant of the overlap of the two ground states associated with 
each Hamiltonian and it can be calculated explicitly once all the negative eigenvectors of 
both Hamiltonians are known. For a finite 2n dimensional lattice the Hamiltonians are 
finite size matrices of size ~ V x V where V is the 2n dimensional volume and their 
eigenvectors can be readily calculated. The resulting chiral determinant has the correct 
magnitude and a phase that exhibits the correct gauge dependence for "smooth" gauge 
fields. For "rough" gauge fileds the phase exhibits a mild breaking of gauge symmetry even 
for anomaly free theories. This problem has been resolved in 0. 

These methods can also be used to formulate a vector theory. In this case the boundary 
conditions along the 2n + 1 dimension are set to be periodic. A Dirac fermion emerges 
with the positive chirality component bound on the wall and the negative chirality bound 
on the anti-wall. If the 2n + 1 dimension is taken to be infinite the two chiralities are 
decoupled and the resulting theory has intact chiral symmetries! Again, this infinity can 
be dealt with the Overlap formalism and now there are no issues associated with the 
phase of the determinant since, for a vector theory, the determinant is real. Therefore, the 
Overlap formalism provides an ideal lattice regularization of vector theories where the chiral 
symmetries are left intact even for finite lattice spacing. Also, the anomalous breaking of 
the axial symmetry is reproduced in an elegant way along with a formula for the index of 
the chiral Dirac operator 0. 

The Overlap formalism was used in a dynamical numerical simulation of the massless 
single and multiflavor Schwinger model with good results ||. Numerical simulations of 
QCD using the overlap formalism would clearly be very appealing. However, as mentioned 
above, such a simulation would require the calculation of all negative eigenvectors of ma- 
trices of size ~ V x V. This makes such a calculation prohibitive for present generation 
supercomputers. 

An obvious alternative (see || and references therein) is to keep the 2n + 1 dimension 
finite and use standard Hybrid Monte Carlo type algorithms to simulate the theory in 
2n + 1 dimensions. Of course, the exact chiral symmetry will be spoiled but it will be 
recovered as the size L s of the 2n + l dimension is sent to infinity. Therefore, even at finite 
lattice spacing one can control the restoration of the regularization induced chiral symmetry 
breaking by using the parameter L s . This involves no fine tuning and, furthermore, since 
the two chiralities decay exponentially away from the wall (anti-wall), one would expect 
that the restoration of chiral symmetry would be exponential i.e. ~ e~ cL % < c. The 
computer cost of such a simulation would be L s times larger than a simulation of standard 
Wilson fermions with the same physical masses. Since for present day supercomputers a 
value of L s greater than 10 — 20 will make simulations impractical, an important question 
to ask is what is the rate "c" of restoration of chiral symmetry and how does it depend 
on the other parameters of the theory and in particular on the lattice spacing. In ||, ||TU| 
some of the issues relating to this question were investigated analytically. Numerical work 
in JII3], [T3| has yielded encouraging results and in particular the interesting work 
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of [13 1 indicates that DWF can successfully address problems related to the evaluation of 
weak matrix elements. However, these works have only marginally addressed this particular 
question. Before full scale dynamical QCD simulations are performed this question should 
be answered. In this paper this question is investigated in the context of the two flavor 
lattice Schwinger model. 

A useful variation of the wall, anti-wall model studied in [14] was proposed in the context 
of vector lattice gauge theory in ||, jpD| . There, instead of having a mass that changes 



sign in two places along the 2n + 1 dimension (say at and L s /2), the mass is kept fixed 
to some positive value mo, but the boundary conditions are taken to be free at the ends of 
the 2n + l dimension. Again, two zero modes with opposite chiralities emerge, but they are 
now bound at the opposite ends of the 2n + 1 dimension and are separated by a distance 
L s rather than L s /2 as in the original model. Therefore, the expectation is that for the 
same L s this model will achieve better restoration of chiral symmetry. 

Another feature added in this model is the introduction of an explicit chiral symmetry 
breaking term that connects the two ends with strength rrif. This gives mass to the fermion 
in addition to the one resulting because of the finite extend L s . The reason for adding this 
term is that it provides linear control over the fermion mass instead of the exponential one 
provided by L s . Furthermore, in a numerical computation, it makes much more sense to 
vary rrif rather than L s in order to control the mass. Therefore, for a given rrif one would 
like to keep L s large enough so that it does not affect the fermion mass in any significant 
way. This method and the associated Overlap formalism will be used throughout this 
paper. 

The theory has five parameters. The first two are the lattice spacing a and the physical 
extent I along one direction of the 2n dimensional box (they are controlled by the pure 
gauge coupling g and the size in lattice units L). The remaining three parameters m , L s 
and rrif all control, to some extent, the amount of chiral symmetry breaking and therefore 
the effective fermion mass. For L s — > oo the theory is chirally symmetric except for the 
explicit breaking introduced by rrif. As a result, the effective fermion mass vanishes linearly 
with vanishing rrif [10|. But for finite L s this is not the case. As mentioned above, even 
for rrif = the restoration of chiral symmetry is expected to be exponential ~ e~ cLs . 

One would expect that the exact continuum solution of the Schwinger model would be 
useful to compare with results obtained on the lattice. Unfortunately, this is only partially 
true. The regularization with DWF introduces to the two dimensional action a four-Fermi 
term with some coefficient. Since, for the two dimensional model, this operator is not 
irrelevant, the continuum theory will be different from a continuum theory with no four- 
Fermi term. Although the continuum theory has been solved with such a term present | J15[| , 
the results can not be directly compared since the value of the coefficient arising from the 
DWF has not been calculated. This problem was encountered in || and [[/J and made the 
comparison with continuum results complicated. Fortunately, for the purposes of this work, 
the continuum results are not needed. In fact, there is a much more relevant comparison 
that can be made. At every step the value of any observable at finite L s can be directly 
compared with its value at infinite L s calculated using the Overlap on the same lattice size 
and lattice spacing. 

The paper is organized as follows: In section |2| the model and the corresponding Overlap 
implementation is reviewed. In section |3] the definitions of the various observables used 
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in this paper are given. In section |] the free theory for finite L s is discussed, the full 
expression for the propagator is given and the "effective" bare fermion mass is identified. 
In section [| some general considerations regarding the interacting theory are presented. 
These considerations lead to specific predictions. In section |^ the results of a full dynamical 
simulation using the Overlap with non zero mass are presented. These results, interesting in 
their own right, are used to compare with the finite L s results of the next section. Section [7| 
describes the results of a dynamical simulation of the 2 + 1 dimensional system for various 
values of the parameters on a fixed physical volume. The algorithm used is a standard 
Hybrid Monte Carlo (HMC) algorithm. The numerical results confirm the predictions 
made in section |5] and together outline the mechanisms of chiral symmetry restoration in 
the model. Section || contains a summary and conclusions. 



2 The Model 

In this section the model and the corresponding Overlap formalism |J implementation is 
reviewed for the benefit of the reader and in order to establish notation ||, [|10|| . The 
following is for a single flavor. The generalization to more flavors is straightforward. 
The partition function of the single flavor 2n + 1 dimensional model is: 

[dU] f [d&dV] [ [d&d$}e- s (1) 



Un(x) is the gauge field, s) is the fermion field and s) is a bosonic Pauli Villars 
(PV) type field, x is a coordinate in the 2n dimensional space-time box with extent L 
along each of the directions, /z = 1, 2, . . . 2n, and s = 0, 1, . . . , L s — 1, where L s is the size 
of the In + 1 direction and is taken to be an even number. The action S is given by: 

S = 3(g , L, L s , m 0} m f ) = S G (U) + S F (V, , U) + S PV (<S>\ $, U) (2) 

where: 

s G = \Y, ReTr l I - u p] (3) 

is the standard plaquette action with g the lattice gauge coupling. In this paper the 
coupling g is exchanged for the parameter: 

(/U) = (4) 



where I is the physical size of the In dimensional box along one of its dimensions and [i is 
a mass related to the photon mass with: 



m 7=V% ( 5 ) 

where Nf is the number of flavors. With these choices is the physical box size in units 
of /i and jil/L = fta is the lattice spacing in units of //. 
The fermion action is: 

S F = - y(x,s)D F (x } s;x',s')ty(x',s) (6) 
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with the fermion matrix given by: 

D F (x, s; x', s') = 5(s - s')fl(x, x') + p-(s, s')5(x - x') (7) 



lp{x,x') = W[^ + l,)U^x)5{x + fi-x') + {l-^)Ul{x')5{x' + fi-x) 

+ (m - 2n)5(x - x') (8) 



f P R 5{\ - s') - m f P L 5(L s - 1 - s') - 5(0 - s') s = 

^(s, s ') = I P R S(s + 1 - s') + P L 5{s - 1 - s') - 5{s - s') 0<s<L s -l (9) 

[ -m f P R 5{0 - s') + P L 5(L S - 2 - s') - 5(L S - 1 - s') s = L s - 1 

iV = ^ do) 

where m is a 2n + 1 dimensional mass representing the "height" of the Domain Wall. In 
order for the doubler species to be removed one must set < m < 2 0]. However, this 
range is further restricted by the requirement that the transfer matrix along the 2n + 1 
direction be positive 0: 

0<m <l (11) 

The gamma matrices are taken in the chiral basis and are the same as in the last reference 
in H. In two dimensions they are: 



71= , 7,= " , 75= n " (12) 






The PV action is designed to cancel the contribution of the heavy fermions in the large L s 
limit. This is necessary because the number of heavy fermions is ~ L s and at the L s —>■ oo 
limit they produce bulk type infinities 0. There is some flexibility in the definition of the 
PV action since different actions could have the same L s —>■ oo limit. However, the choice 
of the PV action may affect the approach to the L s — > oo limit. A slightly different action 
than the one used in |TUJ is used here. This action is easier to implement numerically and 



for finite L s it projects the ground state of the transfer matrix T better; the projector is 
T Ls instead of T Ls/2 (see below). Also, even for finite L s , it exactly cancels out the fermion 
action when rrif = 1 resulting into a pure gauge theory. The PV action is: 

S PV = J2 &fas)D F [mf = l](x,s',x , ,s , )&(x , ,s') (13) 

x,x',s,s' 



The transfer matrix along the 2n + 1 direction for this model is: 

rri „— a' Ha 



(14) 



where a\a are creation and annihilation operators that obey canonical anticommutation 
relations and span a Fock space with vacuum state |0 >. These operators live on the 
sites of the 2n dimensional lattice and carry spin, color and flavor indices. The left/right 
component decomposition of a is: 

I c \ 



(15) 
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The single particle Hamiltonian H is defined by: 



—H 



B- 1 B~ l C 
C^B- 1 C ] B- l C + B 



(16) 



^ 2n 

B(x, y) = - ]T [2 - U,(x)5(x + jl-y)- + fi - x)] + (1 - m )5(x - y) (17) 



(18) 



with G\ = 1, 02 = i in two dimensions. 

The fermionic and PV effective actions can be expressed in terms of the transfer matrix 

as: 

(19) 



e -s eS [L s ] = det(D F [m / ]) = det(B) Ls Tr \T Ls O{m f ) 
e -sg[L.\ = ( det ( J D F [ m/ = l]))- 1 = (det(B) L °Tr [t l *]) 



(20) 



where the operator 0{mf) implements the boundary conditions and contains all the rrif 
dependence: 

°( m f) = U(cnci + m f c ] n c n ){d n di + m f d ] n d n ) (21) 

n 

For rrif — 1 it is the identity operator and for m/ = is a projection operator to a state 
|0'>: 

|0>=n^l 0/> ( 22 ) 

n 

In the infinite L s limit T Ls becomes a projection operator to the ground state of —a) Ha, 



lim T Ls 

L s — >oo 



-Ao 



>< H | 



(23) 



where |0# > and A are the ground state eigenvector and eigenvalue of d^Hd obtained by 
filling all negative energy states. To get an explicit relation between \0 H > and |0 > let 
R be the eigenvector matrix of the single particle Hamiltonian H . The matrices H and R 
have size N x N where N = spin x color x flavor x V. R can be put in the form: 



R 



P- P+ 
Q~ Q + 



(24) 



where the rows labeled P correspond to the left chirality states (with creation/annihilation 
operators c\c) and the rows labeled Q correspond to the right chirality states (with cre- 
ation/annihilation operators a 1 , d). The ± splitting of the columns corresponds to eigenvec- 
tors with positive /negative eigenvalues. With N ± denoting the number of positive /negative 
eigenvalues the size of the P , Q~ matrices is N/2x N~ and the size of the P + , Q + matrices 
is N/2 x iV + . Then it can be shown that: 



I Oh >=U(Zl P i7,i + d kQ7j\U> 



(25) 



i=i 
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From equations IS , 20 and ^3] the effective action for the fermion and PV fields in the 
L s — > oo limit is given by the Overlap formula: 



-8 e g[L s =oo,m f ] 



-SF s [L a =oo,m f ]-Sr/[L s =oo] 



-< H \O(m f )\Q H > 



(26) 



For rrif = 1, < 0#|(9(1)|0// >= 1, corresponding to a theory with no fermions. The rrif = 
case corresponds to massless fermions and the Overlap takes the special form: 

i2 



-S eB [L s =co,m f =0] _ ^ Q^lQ' > | 



It can be shown that: 



< OfflO' >l 



det(Q" 



(27) 



(21 



If N~ = N/2, Q~ is a square matrix and the Overlap will in general be non zero. However, 
if N~ ^ N/2 then Q~ is not a square matrix and its determinant is identically zero. From 
equations ^ |25| and one can see that this arises because of a mismatch in the filling 
levels of | Oh > and |0' >. In order to obtain a non zero overlap one would need to insert 
the appropriate number of creation and annihilation operators to balance the filling levels. 
In fact these operators are the t'Hooft vertices constructed with lattice fields. Then an 
elegant definition of the topological charge q as seen by the fermions arises, H|: 



q = N~ — N/2 



(29) 



where q is naturally integer valued. 

When rrif ^ use of equations ^T], ^5] and [26] yield explicit expressions for the Overlap 
as a determinant of a matrix that is constructed out of P~ and Q~ . These expressions are 
used in the numerical simulation of the Overlap (L s — > oo) in section ^. For more details 
the reader is referred to M and M, |IU . 



3 Observables 

In this section the definitions of the observables that are measured in this paper are given. 
The operators involved are as in 0, |]I0"|| . 

The 2n dimensional fermion operators of the 2n + 1 dimensional theory are constructed 
out of the 2n dimensional fermion fields if), if) as in 



iP{x) = P R V{x,0) + P L *{x,L s -l) 

^(x) = V(x,L s -l)P R + y(x,0)P L (30) 

In the L s — > oo limit of the theory these operators exactly correspond to insertions in the 
Overlap of the creation and annihilation operators discussed in section 0. This will allow 
explicit comparisons to be made between measurements involving if>, if) in the 2n + 1 theory 
with finite L s and measurements with the Overlap involving the corresponding creation 
and annihilation operators. 

The following is a list of definitions of the observables for any L s and the corresponding 
Overlap expressions. The definitions of the actions are as in section [| but for two flavors. 



Use is made of the fact 10 



det(D F ) = det(D{.) (31) 
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The fermion effective action of the 2n + 1 dimensional theory in a background gauge 
field is: 

e -s^[L s , mf ] = J [dy^e-SF = det (p]r[L 8 ,m f ]D F [L a ,m f ]) (32) 
The PV effective action in a background gauge field is: 



e -sf/[L s ] = J ^d9]e~ Spv = det (D F [L s ,m f = l]D F [L s ,m f = 1 
The fermion effective action in a background gauge field is: 



(33) 



det (D F [L S , 


m f ]D F [L s ,m f }) 


det (p F [L 8 ,mf 


= l]D F [L 8 ,m f = l]) 



S eS [L s ,m f ] = e -S^[L s , mf ]-SP/[L s ] 



e -S cH [Ls=oo, mf ] = < H |O( m/ )|0 H > 

The chiral condensate operator is: 

1 2 — ■ 
•W> = -^77 EE [^r( x Wl( x ) + ^l( x Wr 



x) 



x i=l 



(34) 
(35) 

(36) 



The following observable is related to the chiral condensate in a background gauge field: 
P 1 [L S , m f ] = -1 ]T [D F \x, L s - 1, 2; x, 0, 2) + D F 1 (x, 0, 1; x, L s - 1, 1) 



(37) 



Pi[L a = oo,m f ] = ~E[< 0^|c^C(m / )c :r |0 Jf >< 0*10(771,) |0 H > +(c -> d)] (38) 



The t' Hooft vertex operator is: 



to 



= -E 



n (^wi(^)) + n 



X 



(39) 



£ U=l 1=1 

The following observable is related to the t' Hooft vertex in a background gauge field: 



P 2 [L s ,m f ] = J [d^d^] w e 
P 2 [L s ,m f ] = 1 [D F 2 (x, L s -1,2;x, 0, 2) + D F 2 (x, 0, 1; x, L s - 1, 1) 



P 2 [L S = oo, m,] = i £ [< 0* ^(m,)^ > 2 +(c -> d) 



(40) 



(41) 



The observables e Scff , Pi and P 2 are interesting because they correspond to Overlap expres- 
sions and, therefore, their values at the L s = oo limit are calculable. Furthermore, these 
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quantities are sensitive to the topology of the background gauge field. The expectation 
value of the fermion condensate is: 



< ^ >= / [dU] I [dVdV] 



< ipip > 



ipip e s 

fdU±E x < H \clO{m f )c x \0 H >< H \O(m f )\0 H > +(c 



(42) 



d) 



L s =oo~ 



So 



JdU < H \O(m f )\0 H > 2 e~ s o 
The expectation value of the t' Hooft vertex is: 



< w >= I [dU] I [dtydV] 



w e 



(43) 
(44) 

JdU < H \O{m f )\0 H > 2 e~ s ° ^ 

Notice that a numerical evaluation of < ipip >l s=00 and < w >L a=00 requires two separate 
pure gauge simulations, one for the numerator and one for the denominator, i.e. the fermion 
determinant is treated as an observable ||, f7|. 



< w > Ls 



I dU± Ex [< 0^140(7^)410^ > 2 +(c -> d) 



4 The free theory 

In this section the propagator is given, its singular part is identified and the bare fermion 
mass which is a function of L s ,m ,mf, is extracted. As a verification of this result the 
smallest eigenvalue of the free In + 1 dimensional Dirac operator is also calculated. 

The propagator has been calculated for the infinite L s case in the first reference of 0. 
The propagator has also been calculated for the model described in section || but only in 
the limit where exponentially small contributions in L s could be ignored |J. The size of 
these contributions was alluded to but no explicit expression was given. Since in this paper 
the interest is on the behavior of these contributions, the full calculation is worked out. 
Because the general form of the propagator is the same as in 0, an effort has been 
made to keep similar notation. 

The free 2n + 1 dimensional Dirac operator of eq. ^ in momentum space is: 

D F (p : s, s') = vp5(s — s) — b(p)8(s - 

and 

2n 

hip) = E [1 

p M = sin(p M ), n=[l,...,2n] 
M(s, s') =8(s + l- s') - 



s') + ^^M{s, s') + - — -M^ (s, s') (46) 



- cosQ^)] + 1 - m 



(47) 



2iik 



£ k, = [0,1,. ..L-l] 



-S(s' -0)6(L a -l-s)(l + m f ) 



(48) 
(49) 
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where the 5 functions are understood as having "period" L s . Notice that the first two terms 
of eq. |46| are the same as for Wilson fermions with mass (1 — m ), where < m < 1 (see 
sect. 0). The second order free Dirac operator is diagonal in spin: 



(50) 



where Q± have no spin indices and: 



Q+{p : s, s') = [-6(p) + M(s, s')} [-b{p) + M\s, s')} + f 
VL-{p : s, s') = fi+(p : L s — 1 — s, L s — 1 — s') 



(51) 



where: 



2n 

—2 \ —2 
p - \ .. 



(52) 



The inverse of the second order free Dirac operator DjDp must therefore be of the 
form: 



G = l -±^G + + l^G. 



D f D 



G 



2^2 

where G± have no spin indices. From general considerations G± must be of the form 
G + {p : s, s') = A e- als ~ s ' 1 + A^ 3 ^ + ^ 2e — 



(53) 



+ A 



e -a(L s -l+s-V) _|_ e -a(L 8 -l+s'-s) 



G-(p:s,s') 



G + {p : L s - 1 - s,L s - 1 



(54) 
(55) 



The coefficients in eq. |55| are momentum dependent. The coefficient "a" is the solution 
of the equation: 

1 _|_ 52 _^ p2 

cosh(a) = (56) 

2b 

A straight forward calculation results to the following expressions for the remaining coeffi- 
cients: 

A 



where, 



1 



26sinh(a) 



A x = -^{b-e~ a ){l -mj) 



An 



B 
'A 
B 



(e a -b){\ -m)) 



A, 



B 



2bm f sinh(a) + e - a{Ls ~ 1] (e- 2a [e a - b] 



m. 



,2a 



2 f„a 



e a — b) + [4mj6sinh(a)] e 



-a(L s -l) 



+ 



m 2 f {b-e~ a )+e 



2a / a 



6) 



-2o(i a -l) 



(57) 

(58) 
(59) 
(60) 



(61) 
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In order to identify the singular part of the propagator an expansion in the variables p, 



rrif and (1 — m 



.(£.-1) 



is done treating these variables as numbers with magnitudes much 



smaller than one. The only singular amplitudes are A 2 and A m and the resulting expression 
for the singular part of G± is: 

1 



G 



singular 



(p : s, s) 



p 2 + m^ ff 



{m (2 - m )e 
- m eff (l - m 



a(L s — 1— s+L s — 1— s') 



-a(L s 



-1+s— s') _|_ g-a(L s 



-l+s'-s) 



with, 



ra cr r = rn (2 - m ) m/ + (1 - m ) J 



} (62) 



(63) 



The expression for Q s ^ ngular can b e obtained from eq. B2| and [55|. The propagator for the 



free theory is then given by: 

D F l = D F G (64) 

To verify the above result the smallest eigenvalue of D F D F is also calculated. The 
leading order term in an expansion as the one used above is: 



A, 



p 2 + m 2 eS + 0(4) 



(65) 



and is proportional to the inverse of the singular part of G± as it should be. 

The interesting result of the above analysis is eq. This is the mass of the lightest 
mode of the theory and is controlled by rrif, mo, L s . This formula strongly suggests the 
pattern for chiral symmetry restoration in the model. 



5 The Interacting Theory, General Considerations 

The question on how the chiral symmetry restoration rate depends on the various parame- 
ters of the model ultimately can only be answered when the full theory is considered. This 
is done using numerical simulations in section [?]. However, some general considerations that 
hint to the expected behavior are useful in order to guide the numerical experiments and to 
provide an understanding of the results. Such considerations leading to specific predictions 
are presented in this section. These predictions are confirmed by the numerical simulations 
and therefore sketch the mechanisms of chiral symmetry restoration in the model. 



5.1 The Effects of Topology 

In this section a special kind of trajectory [17|, || that "cuts" through the various topolog- 
ical sectors of the configuration space is considered. The values of the relevant observables 
of section |3| are calculated along this trajectory. 

Consider the following set of U(l) gauge link U configurations labeled by the continuous 
parameter r 6 5ft |fTT| , 0: 



Ui(ni,n 2 ) 
U 2 (n 1 ,n 2 ) 




if rii L 
if rii = L 



exp 



(66) 
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where Ux, U2 are links defined on the sites nx, ^2 of the two dimensional torus of size Lx L. 
This configuration is periodic along the second direction but has a discontinuity along the 
first. When r is an integer the electric field strength is uniform E = 2ttt/L 2 . When r is 
not an integer E has a discontinuity at (nx — 1, n 2 — 1). The topological charge of the gauge 
field is defined as ||T6| : 

where the sum is over all plaquette variables U p and the logarithm has log(l) = with the 
cut along the negative real axis. It is straightforward to verify that q is an integer that 
changes values as r is varied between integer values (see figure 1). For more information 
on these configurations the reader is referred to [[H], ||. Here it is worthwhile to point 
out that this trajectory in configuration space is interesting because it connects the trivial 
configuration with uniform configurations that have non zero topological charge. As such, 
the configurations at integer r can be thought of as local minima and therefore in some 
sense they represent vacua of different topological charge. The path that connects the 
configurations of integer r is certainly not unique but it can nevertheless provide insightful 
information on how the transition between sectors of different topological charge takes place. 
Therefore, although one can not extract from this trajectory quantitative information, the 
qualitative information will be quite useful if, in particular, similar features are observed 
when the full configuration space is considered in the HMC simulation of the model in 
section |7|. 

In the following, the observables in equations [32] - 145| are considered in the presence of 
the gauge field configuration of eq. The boundary conditions for both fermion flavors 
in the two dimensional space are taken to be antiperiodic. Of course, for the full dynamical 
U(l) theory, the choice of boundary conditions is irrelevant for as long as they are the same 
for all flavors (they can change by U(l) phases at the boundary). 

The L s = 00 quantities are calculated using the Overlap formula. The finite L s quan- 
tities are obtained by explicit computation of the determinants and inverses of Dp. The 
reader is reminded that e - s c«[ L s=^,m f =o\ j g no ^. zero on jy j n ^ e topological sector q = 0. 

The quantity related to the chiral condensate P\[L S = 00, rrif = 0] is identically zero in all 
sectors. This is a consequence of the exact SU(2)l x SU(2)r chiral symmetry. Therefore 
a non zero Px[L s ,mf = 0] at some finite L s will indicate explicit breaking of the chiral 
symmetry. The quantity P2[L S = 00, rrif = 0] related to the t' Hooft vertex is not zero only 
in sectors q = ±1. This represents the anomalous breaking of the U(l) axial symmetry. 

In figure 1 the topological charge along with Sq, S^, , S e g and S e ff + Sq is shown 
for L = 6, fil — 3.0, m = 0.9, rrif = 0, and L s = 14 along the trajectory of eq. |66|. The 
absolute scale in these figures is of course irrelevant. One can already see the separation 
of the different topological sectors with the q = sector as the absolute minimum and the 
other sectors as relative minima. It is interesting to observe how the addition of and 
Sq changes the fermion action S^ s . 

In figure 2 S e s, P\ and P2 are shown along the same trajectory and for the same 
parameters but for various values of L s . In the S e g figure at r = 2 the curves with larger 
S'eff correspond to larger L s with L s = 4, 6, 8, 10, 12, 14. The curve that sharply increases 
to 00 at r ~ 0.5 corresponds to L s = 00. It is clear from this figure that in the q = sector 
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the L s = oo limit is approached very fast. The q 7^ sectors carry large actions relative to 
the q = action and therefore they also approximate the L s = oo limit very well. The only 
regions that suffer from a slow approach to the L s = oo limit are the two regions in the 
immediate neighborhood where r changes from to ±1. There, the L s = oo result raises 
discontinuously but the finite L s results approach this discontinuous behavior slowly. 

In the Pi figure at r = the curves with smaller P\ correspond to larger L s with 
L s = 4, 6, 8, 10, 12, 14. The L s = oo curve is identically zero for all r and is not plotted. 
From this figure it is seen that Pi in the non zero topological sectors is negligible even at 
L s = 4. In the zero topological sector it decreases very fast with increasing L s for all values 
of t except for the same two regions where r changes from to ±1. There, two large spikes 
appear with heights that essentially do not change with increasing L s . Instead, their width 
slowly shrinks with increasing L s . As a result, chiral symmetry is slowly restored in these 
regions. 

In the P2 figure at r = 1, the curves with larger P 2 correspond to larger L s with 
L s = 4, 6, 8, 10, 12, 14, 00. The L s = 00 curve can be distinguished by the discontinuous 
behavior at the places where q changes form to ±1 and from ±1 to ±2. From this figure 
it is seen that P2 approaches the L s = 00 limit very fast for all r except again for the 
regions where q changes. 

The slow approach of S e ff and Pi to the L s = 00 limit in the regions of changing 
topological charge is not particularly troubling since the effect is a small percent of the 
values that they acquire along the trajectory. Not so for Pi. The two "spikes" present a 
large size contribution to a quantity that should otherwise be identically zero. In figure 
3 Pi is plotted vs. L s for various "cross sections" of figure 2. For r = (diamonds) Pi 
decays exponentially. For r = 0.3, 0.4 (squares, crosses) the rate of the exponential decay 
decreases. Finally, close to the point where the configuration changes topological charge, 
r « 0.55 (octagons), there is almost no decay at all. 

The slow approach to the L s = 00 limit in the neighborhoods of changing topological 



charge q is expected. According to the topological charge definition eq. 29, q changes 



because the number of negative eigenvalues of the single particle Hamiltonian H of eq. 16 



changes. Therefore, in that neighborhood there are configurations for which H has zero 



eigenvalues. These configurations were identified in || and also discussed in For 
these configurations the Overlap formulation is not well defined because the ground state is 
degenerate and, in turn, the finite L s theory experiences large correlations along the 2n+l 
dimension. However, the set of configurations for which H has an exact zero eigenvalue is 
of measure zero. As a result the Overlap is a well defined formulation. On the other hand 
regions of the configuration space surrounding these special configurations are characterized 
by small decay rates. The importance of these regions is determined dynamically and 
only a full simulation of the theory can accurately probe their effect to the exponential 
decay. However, close to the continuum limit one would expect that these regions are 
severely suppressed since the pure gauge action separates sectors with topological charge 
that differs by one unit with barriers of energy ~ (see Sq in figure 1). But away from 
the continuum limit, where most numerical simulations are performed, these regions may 
become important and contribute to the explicit breaking of chiral symmetry. 

In order to gain some understanding about this mode of chiral symmetry breaking the 
expectation value of the chiral condensate < ipip >, eq. 03L was calculated in a configuration 
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space restricted only on the trajectory of eq. i.e. < ipip > was calculated as in eq. ^ but 
with the path integral over the gauge field configuration space replaced by an integral over 
the single trajectory of eq. ^ parametrized by r. This calculation was done numerically. 
As can be seen from figure 2, the integrands are relatively smooth and become negligible as 
\t\ is increased above 2. Therefore, it was enough to choose a reasonably fine discretization 
and sum over a finite range only. The summation was done for step size dr = 0.025 and 
over the range \t\ < 2.5. < ipip > is plotted vs. L s in figure 4 for an L = 6 lattice, m = 0.9, 
rrif = and various gauge couplings fil/L, fil = 3.0,2.5,2.0,1.5,1.0 (diamonds, squares, 
crosses, octagons, stars). For fil = 3.0 there is decay with a fast rate until L s = 10. For 
L s = 12 and above there still is decay but with a smaller rate. As can be seen from this 
figure the functional form of the fast decay is exponential. The slower decay is consistent 
with exponential but it turns out that it is also consistent with power law or exponential 
times power law behavior. Since the functional forms are further "contaminated" by higher 
order effects in order to be able to distinguish between the different types of decay it is 
estimated that calculations up to L s w 30 will be needed. This is beyond the purpose of this 
simple test and the available computer resources. In any case, the important observation 
is that as the gauge coupling is decreased the inflection point moves to larger L s and the 
slower of the two decays becomes faster until at \xl = 1.0 there is no visible inflection below 
L s = 18. This phenomenon can be easily understood by looking at Pi in figure 2 (/z/ = 3.0 
there). For L s smaller than about 10 the explicit breaking that occurs inside the q = 
sector dominates and the contribution of the "spikes" to the expectation value is small by 
comparison. When this breaking has almost completely disappeared, L s > 10, the breaking 
that comes from the "spikes" dominates. This breaking disappears in a slower fashion as 
the width of the "spikes" shrinks. However, as the gauge coupling is decreased the regions 
of changing topology where the "spikes" are located are weighted less by the pure gauge 
action and the slower decay is overshadowed by the initial faster exponential decay. 

Therefore, one sees that there are two distinct mechanisms that control the restoration 
rate of chiral symmetry. One is related to the restoration in the q = topological sector. 
The other is related to the topology changing regions of the gauge field configuration space 
and in particular to the regions that connect the q = and q = ±1 sectors. 

To conclude this section the effects of m/ are presented. The observable Pi along the 
trajectory of eq. ^6| is plotted in figure 5 for L = 6, \il = 3.0, m = 0.9, rrif = 0.1 and 
L s = 4, 6, 8, 10, 12, 14, 16, oo. At r = the curves with larger Pi correspond to larger L s . It 
is clear that the L s = oo limit is approached rapidly and with no complications. Again the 
region of changing topology approaches the L s — > oo limit slowly but it is away from it only 
by a small percentage. One would expect that the non zero rrif behavior will persist down 
to some small value of rrif before signs of the rrif = behavior described above appear. 
One can visualize how figure 5 changes with decreasing rrif. As rrif becomes smaller Pi 
tends to zero in the various topological sectors. However, its value, at the place where the 
topological charge changes from to ±1, remains roughly constant resulting in the spikes 
of figure 2 . 

In order to see more clearly how the L s = oo limit is approached for rrif = 0.1, the 
expectation value of the chiral condensate < ipip >, eq. |4^, is calculated in a configuration 
space restricted only on the trajectory of eq. for various L s and for L s = oo. In figure 
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6 < ipip >l s — < ip4> >oo / < >L a is plotted vs. L s . Again, one can see that there 
is an inflection at about L s = 10. However, when the inflection occurs the L s = oo value 
has already being approached to better than 0.3%. This behavior is encouraging as far as 
numerical simulations are concerned. Since most numerical simulations are done for small 
but non zero masses one would expect that for some range of masses the effects of the 
topology changing configurations to these simulations will generally be small and perhaps 
even lost in the statistical noise. 



5.2 The Effects of Gauge Field Fluctuations 

In this section the effects of gauge field fluctuations to the chiral symmetry restoration rate 
are discussed. 

The two mechanisms of chiral symmetry restoration identified in section |5.1| will be 
affected when the gauge field is allowed to fluctuate. The effect of dynamical gauge fields 
to the mechanism that restores chiral symmetry in the zero topological sector can be seen 
by measuring < ipip > in an ensemble of configurations generated by applying small fluc- 
tuations to the trivial configuration. In particular, consider gauge field configurations with 
links U = e ir7T where r is a random number in the range — e < r < e with e a small number 
that controls the size of the fluctuations. These configurations have a "flat" distribution 
and in this ensemble < ipip > is obtained by calculating < Pi > / < e~ ScB >. In figure 7 
< ipip > is plotted vs. L s for various values of e = 0.4, 0.3, 0.2, 0.1 and 0.01 corresponding to 
octagons, stars, squares, crosses, and diamonds. The value at each point was calculated in 
an ensemble consisting of 40 configurations. Antiperiodic boundary conditions have been 
used for the fermions and L = 6, rrif = 0, m = 0.9. It can be seen that as the size of the 
fluctuations decreases the chiral symmetry restoration rate increases. 

Some insight to this behavior can be gained by considering the following "heuristic" 
argument. The Dirac operator of equation [?| can be rewritten as: 

D F = ^naivc + M 



pa a i ve (x,x') = \ \Un(x)6(x + fi-x) - Ul(x')5(x' + p, 



X 



M = -B + ^~y^M + — (68) 

where B is the Hermitian matrix given in eq. [T7] and M is given by equation For 
the free theory in momentum space D F = vp + Ai(p) where Ai(p) can be read from 
eq. One can think of A4(p) as being a momentum dependent mass matrix. The 

smallest eigenvalue of AiAi^ is obtained at zero momentum and is equal to with 
m e ff = m Q (2 — mo) m/ + (1 — m ) is , eq. |63]. The quantity (1 — m ) is the smallest 
eigenvalue of the matrix B. When interactions are turned on the smallest eigenvalue of B 
will shift to values (1 — m' ) larger than the free theory (1 — mo) value for the simple reason 
that the matrices U^(x)S(x + fi — x') that make up B are unitary. Therefore, the smallest 
eigenvalue of AiA4^ will now be m' eS 2 , m' efi = m (2 — m ) rrif + (1 — m' ) Ls . For small 
fluctuations one may still be able to think of M. as a mass matrix. If that is the case the 
lightest mass in the theory will be m' e{i and one sees that for rrif = the chiral symmetry 
restoration rate will become faster as the fluctuations become smaller. 
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Finally, the quantitative effect of dynamical gauge fields to the mechanism related to 
topology changing is complicated since it involves some understanding about the volume of 
configuration space that contains gauge field configurations for which the Overlap Hamil- 
tonian has near zero eigenvalues. As already mentioned, one would expect that this volume 
will shrink as the continuum limit is approached because the pure gauge action introduces 
energy barriers of size ~ l/^o between sectors with topological charge that differs by one 
unit. 



5.3 The range of m 

In this section the allowed range of mo and the effects it has in the approach to the 
continuum limit are discussed. 

In order for the doubler species to acquire masses of the order of the cutoff m must 
be in the range < m < 2 [|Tj. As it was found in |5]], in order for the transfer matrix 
to be positive definite for all gauge fields, the range of mo should be further restricted to 
< itlq < 1, eq. |TT| For QCD, any value of mo in this range should lead to the same 
continuum limit since the local and global symmetries of the theory remain unchanged. 
On the other hand, different choices of mo, will result in different ways of approaching this 
limit. For the Schwinger model there is an additional complication. As mentioned in the 
introduction, a four-Fermi term is a marginal operator in two dimensions. Although it is 
not explicitly introduced in the action, the DWF regularization introduces such a term 
with a coefficient that depends on m [§], @- Therefore, different values of mo will lead 
to different continuum limits. For this reason, a quantitative study of the effects of m to 
the approach to the continuum limit in the Schwinger model is complicated and will not 
be done here. However, there are some important generic features that can be discussed. 



Consider the Overlap Hamiltonian H, eq. [16] as a function of the 2n + 1 dimensional 
mass mo H. It is easy to see that for mo = — oo, H has the same number of positive 
and negative eigenvalues, N + = N_. It is also easy to see that for m < 0, H can not 
have a zero eigenvalue. As a result H has q = 0, for m < 0. As mo is increased from 
zero, H can develop zero eigenvalues and as a result q ^ (eigenvalues cross zero altering 
N~). Close to the continuum it can be easily seen that most crossings will occur around 



m o = 1. Farther away from the continuum pure gauge numerical simulations [IS show 
a finite region [mo min , mo max ] where most crossings occur. This suggests that one should 
keep m 0max <C m for two reasons. First, one would like to be as far away as possible from 



the region where crossings occur since, as discussed in section |5TT| , it is there that the decay 
rates become small. Second, as can be seen from the definition of the topological charge, 
eq. it is only then that the full effects of non-trivial topology are visible to the fermions. 
In particular, one would expect that by keeping mo max < mo quantities that are sensitive 
to topology will have a "smoother" approach to the continuum limit. 

From the above discussion it appears that at any coupling the safest choice would be 
to set mo to its largest allowed value mo ~ 1. Even then, because in general < mo min < 
mo max < 2 it may be that far away from the continuum 1 < m 0max => m < m 0max . If this 
situation occurs the finite lattice spacing errors at that coupling will be large and the chiral 



symmetry restoration rate will be slow. However, the numerical simulations of [18] indicate 



that even for reasonably strong couplings m 0max < 1. The couplings used in this paper fall 



17 



in this category (see figure 1 in [18]) and mo is kept fixed to tuq = 0.9. 

In ]13| it was found that for QCD some tuning (mo = 1.7) was necessary in order to 
sufficiently restore chiral symmetry at L s = 10. However, the configurations used there 
were not generated with the DWF action but rather with the staggered fermion action at 
j3 = 5.7. It is possible that this large value of mo is at least due in part to the "semi- 
quenched" nature of the calculation and if the configurations are generated with the DWF 
action setting m « 1 may be sufficient. 



6 Dynamical Simulation of the Massive Theory with 
the Overlap 

In this section a full dynamical simulation that measures the chiral condensate < ipip > 
and the t' Hooft vertex < w > of the two flavor massive Schwinger model using the Overlap 
formalism of the model of section [2] will be presented. The result is interesting on its own 
right and it also provides the L s = oo numbers that will be used to compare with the 
results of the dynamical simulation for finite L s presented in the next section. 

The expectation values < ipip > and < w > can be calculated numerically using equa- 



tions |43| and |45] respectively. As it is evident from these equations a pure gauge simulation 
with action Sq must be performed with the overlap factors appearing in the numerator 
and denominator treated as observables. The pure gauge theory expectation values of the 
overlap factors in the numerator and denominator are then divided to produce < ipip > 
or < w >. A heat bath algorithm was used to generate configurations with the standard 
Wilson plaquette pure gauge action of eq. 

Measurements of < w > for the massless, rrif = 0, theory were performed in ||. The 
reader is referred there for more details on the method and results. However, the result 
for w in that reference can not be compared directly with the rrif = result here since a 
different implementation of the Overlap was used. Furthermore, the single plaquette action 
used there was a heat kernel action while the standard Wilson plaquette action is used here. 
The result for < w > /m^ vs. L s on an L = 6, mo = 0.9, lA = 3.0 lattice is given in figure 
8. This quantity acquires a non zero vacuum expectation value even for m/ = (dotted 
lines) as a result of the anomalously broken U(l) axial symmetry. 

The result for < ipip > /m 7 vs. rrif on the same lattice and for the same parameters is 
presented in figure 9. The behavior of < ipip > vs. rrif is interesting in that < ipip >~ rrif 
for rrif < 0.1 while < ipip >~ rrif for rrif > 0.1. Fits to Am P f for rrif < 0.1 and for 
rrif > 0.1 are shown in the same figure. Both fits have a x 2 P er degree of freedom of about 
one. For m f < 0.1 p = 0.996(3) while for m f > 0.1 p = 0.32(2). 



This type of behavior was found by analytical continuum calculations in [19]. In par- 



1/3 

ticular, the linear behavior was found to take place for rrifL 1 while the rrif behavior 
was found to occur for rrifL ^> 1. Unfortunately, the coefficients calculated in that refer- 
ence can not be directly compared with the ones here because, as previously discussed, an 
additional four- Fermi interaction is induced by the DWF regularization ||, ]/J. This term 
in two dimensions is not irrelevant and it will contribute to the continuum limit. For the 



same reasons < w > can not be directly compared with the continuum results of [19[ 
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7 Hybrid Monte Carlo Simulation 



In this section a full dynamical simulation of the two flavor Schwinger model for finite L s 
is presented. The algorithm used is the standard Hybrid Monte Carlo (HMC) algorithm 
|20|| . The expectation value of the chiral condensate < tpip >, eq. ^2], is calculated and 



used to monitor the amount of explicit chiral symmetry breaking. The reason for using 
< iptp > instead of the pion mass is that in two dimensions there is no spontaneous chiral 
symmetry breaking and although the pion mass vanishes for vanishing fermion mass the 
pion is not a Goldstone particle. In this sense < ipip > is as good a probe of chiral symmetry 
breaking as is the pion mass, but unlike the pion mass it has the practical advantage of not 
requiring large lattice size along the time direction in order to measure the decay of the 
pion correlator. The effects of the anomalously broken U(l) axial symmetry are monitored 
by measuring the expectation value of the t' Hooft vertex < w >, eq. §3. 



The smaller size lattices were simulated on the workstations of the Physics Department 
of Columbia University. The larger lattices were simulated on the Silicon Graphics Power 
Challenge Array computer system at NCSA UIUC and also on the C90 supercomputer at 
PSC. 

7.1 The algorithm 

The HMC Hamiltonian is given by: 

#hmc = ^P 2 + S G [U]+x t (D^[m f ,L a ]D F [m f ,L a ])~ 1 x 

+ &(pl[m f = liL s \D F [m f = liL s \)<& (69) 

where P is the HMC momentum, Sq is the pure gauge plaquette action given in eq. ||, U 
is the £7(1) gauge field, x is the pseudofermion field, $ is the bosonic PV field and D F is 
the three dimensional Dirac operator of eq. ^. 

The HMC trajectory length r is set to r = 1, the step size is t/N and the number 
of steps N is adjusted according to the size of the effective bare fermion mass. Typical 
simulations with rrif > 0.05 are done with 50 HMC steps. For the rrif = simulations the 
number of steps is set to 100 — 400 depending on the size of L s . This is necessary because, 
for rrif — 0, trajectories that cross between different topological sectors experience large 
HMC fermion forces (for more details see section [7. 2|) . The Conjugate Gradient residual is 
set to 10 -8 . Typically, the number of Conjugate Gradient iterations is around 50 and does 
not exceed ~ 100 for rrif = and L s = 14. The number of measurements is ~ 8, 000 except 
for the more "expensive" 12 x 12, 10 < L s lattices where 1, 000 — 2, 000 measurements were 
performed. 



7.2 HMC and topology 

In this section problems related to the sampling of non zero topological sectors with the 
HMC algorithm at small fermion masses are discussed. 

The HMC algorithm is very successful provided the fermion mass does not become 
very small. If the fermion mass becomes very small then the effects of topology are not 
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reproduced correctly. In particular, the measurement of < w > becomes problematic. If, 
for some fixed volume, the fermion mass is made very small then similar analysis as in 



21f] indicates that the effect of the zero modes coming from different topological sectors 
becomes important. In particular, as can be seen from the Overlap implementation 0, 
||, at finite physical volume, at zero mass and in the L s —>■ oo limit, the operator w 
receives contributions only from sectors ±1, while the fermion Boltzman weight (fermion 
determinant) is not zero only in sector 0. As mj is turned on (and/or L s is decreased 
from infinity) the fermionic determinant becomes non zero in sectors other than and the 
operator w receives contributions from sectors other than ±1. Because of this behavior at 
small fermion mass the HMC algorithm will mostly sample the sector where the observable 
w receives small contributions. The algorithm will infrequently visit the sectors ±1, but 
when it does the observable w will receive large contributions to "make up" for the small 
sampling rate. As a result, when the fermion mass is decreased a larger number of HMC 
iterations will be needed to sample the ±1 sectors correctly. Therefore, for a fixed amount 
of computer time, if the fermion mass becomes very small the important contributions 
may not even be sampled at all and as a result not only the expectation value of w will 
be underestimated but also the associated error. This type of difficulty has already been 
noticed in simulations of the Schwinger model [22|, \T2\. 

The problem described above leaves a clear signature in the time history of w. In 
figure 10 the time history of w is given for six different fermion masses at L s = 14. As 
the fermion mass is decreased, the average around which w fluctuates decreases. This 
decrease is compensated by the large contributions received from the q = ±1 sectors. These 
contributions start to appear as "spikes" in the time history. The smaller the fermion mass 
the larger the "height" of the spikes but since the corresponding Boltzman weight becomes 
smaller their frequency also decreases. Notice the different scale of the rrif = 0.05 and 
irtf = 0.0 graphs. 

But it is not only < w > that is affected at small fermion masses. If, along an HMC 
trajectory, the topological charge changes then the fermion determinant changes by a large 
amount. As a result, the HMC fermion force becomes large and the HMC step size errors 
become large. Therefore, as the fermion mass is made smaller one must adjust the HMC 
step size accordingly. Finally, it should be stressed that these difficulties are not related to 
DWF but rather to the HMC algorithm. In particular, the simulations using the overlap in 
section |6| do not suffer from these problems since there the fermion determinant is treated 
as an observable in a pure gauge heat bath. 



7.3 m f = 

In this section a strict test of chiral symmetry restoration is done. 

A full dynamical simulation using the HMC algorithm is performed with the explicit 
fermion mass rrif set to zero so that the only breaking of chiral symmetry comes from the 
finite extent L s . The restoration rate at fixed physical volume and various lattice spacings 
is studied by measuring < ipip > for various values of L s . Although this provides a strict 
test one must keep in mind that the non-zero topological sectors may be suppressed more 
than they should for the reasons mentioned in section |7.2| . This means that the rate of 
restoration of chiral symmetry observed here is mainly due to effects that occur in the zero 
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topological sector and its vicinity. The following results were obtained at fixed mo = 0.9. 

In figure 11, < ipip > /m 7 is plotted in a "log" plot vs. L s at fixed physical volume 
fil = 3.0 and for various lattice spacings fil/L = jia where L = 6, 8, 10, 12 corresponding to 
the lines from top to bottom. Data for L = 4, L s = 6 — 10 are statistically indistinguishable 
from the L = 6 data and are not plotted. For L s = 6 — 10 the decay is consistent with 
exponential with a rate that becomes faster as the lattice spacing decreases. For L s = 12—22 
the decay is again consistent with exponential but with a slower rate. Again, this rate 
becomes faster as the lattice spacing decreases. Also, the percent change of the rate at 
L = 6 is ~ 54 ± 6% but at L = 12 is ~ 31 ± 6%. The fits shown are two parameter fits to 
< ipip > /m 1 = Be~ cLs . The x 2 P er degree of freedom is smaller than one for all the fits 
except for the L = 12, L s = 6 — 10 data that have a x 2 P er degree of freedom of ~ 3. 

The exponentiated rate e~ c of the various fits vs. 1/L ~ a is shown in figure 12. The 
diamonds correspond to the L s = 6 — 10 fits while the crosses to the L s = 10 — 22 fits. One 
can see that e~ c is roughly a linear function of 1/L ~ a for the L s = 6 — 10 fits and for 
L = 8, 10, 12. However, more data at smaller lattice spacings are needed before one can be 
confident that scaling has set in and that this is the correct scaling form. 

Although the above fits are all consistent with exponential decay, power law decay of 
the form < ipip > /m 7 = BL S ~ P can be excluded with some confidence only for the fast 
decay, L s = 6 — 10, at the smallest lattice spacing, L = 12. A power law fit to this data has 
X 2 per degree of freedom of ~ 32. More statistics and larger L s will be needed in order to 
clearly establish the type of decay for the other cases. For example, at the largest lattice 
spacing, L = 6, it is estimated that the error bars will have to be reduced from their few 
percent size down by a factor of about ten. Alternatively, the error bars can be kept at 
the few percent level but then it is estimated that L s will have to be made as large as 
pa 30. Both approaches are beyond the purpose of this paper and the available computer 
resources. 

The facts that the decay changes for L s > 10, that the slower decay approaches the 
faster one as the continuum limit is approached and that both decays become faster as 



the continuum limit is approached can all be understood from the analysis in sections |TT 



and |5.2|. Finally, the t' Hooft vertex < w > was also measured but, as expected from the 



discussion in section [7.2| , its value and associated error is underestimated. In particular, 
its value is much lower than the corresponding Overlap value of section || For this reason 
that data is uninteresting and is not presented here. 

7.4 m f ^ 

In this section the mf ^ case is studied. Since typical QCD simulations are done for non 
zero fermion mass the results of this section are of practical interest. 

Dynamical simulations are performed with masses m,f large enough, 0.1 < rrif, so that 
the effects of topological sectors q — 0, ±1 are not miscalculated due to problems associated 
with the HMC algorithm as described in section [7.2|. Both < ipip > and < w > are measured 



and their approach to the L s = oo limit is studied and compared with the L s = oo results 
of section This is done for fixed physical volume, and various lattice spacings. The 
parameter m is kept fixed at 0.9. 
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In figure 13 

volume and lattice spacing, id = 3.0, L = 6. 

expects that < ipip > /m 7 ~ wi e fr- Therefore a fit of < ipip > /m 



< ipip > /m 7 is plotted vs. L s for rrif 

According to the results in section 

3 



0.1, 0, 2, 0.3, 0.5 at fixed physical 

one 



vs. L, is made to 



the form A + Be~ s . All fits have a x P er degree of freedom « 1 — 2. In these figures the 
cross is the coefficient A and the dotted lines are the L s = oo result of figure 9 plus/minus 
the error. One can see that as mj becomes larger the L s = oo result is approached faster. 
This is in accordance with naive expectations born out from the free theory formula for 
m e ff eq. 33. In figure 14 < w > /m 2 is plotted vs. L s for the same parameters as in figure 
13. The fits are again to a form A + Be~ cLs and have \ 2 per degree of freedom k 1-2. 
One can see that the L s = oo result has already been approached at L s = 6. 

The effects of changing the lattice spacing at mj- = 0.2 can be seen in figures 15 and 16. 

r — / I 3 

In figure 15 < ipip > jm 1 is plotted vs. L s at fixed physical volume ill = 3.0 for different 
lattice spacings fil/L = fia, L = 4,6,8, 10. The fits are again to a form A + Be~ cLs and 
have a x 2 per degree of freedom of ~ 1 — 2. One can see a similar behavior as the one in 
section |7.3| . As the lattice spacing is reduced the rate c of the exponential approach to the 
L s = oo result increases. For example, < ipip > a t the larger lattice spacing L = 4 decays 
with c = 0.54(3) but at the smaller lattice spacing L = 10 it decays faster with c = 1.1(1). 
However, one should note that for L = 8 and L = 10 the rate saturates and is essentially 
dictated by the L s = 4, 6 points with the L s = 6 point very close to the L s = oo result. If a 
second slower rate sets in for 10 ^ L s it is unimportant and is lost in the statistical noise. 
Some insight to this behavior can be gained from the analysis at the end of section |5.1| (in 
particular see figure 6). Similar behavior is observed in figure 16 for < w > /ml.. All fits 
have x 2 P er degree of freedom rs 1 — 2. 

Similar results are obtained if one keeps the physical volume and m/ in physical units 
fixed while changing the lattice spacing. This can be seen in figures 17 and 18 by comparing 
the L = 10, rrif = 0.2 data (diamonds) with data at L = 4, rrif = 0.5 (squares). In these 
graphs the physical volume is fixed at fil = 3.0, and rrif in physical units is fixed at 
rrifL = 2.0. Again the decay rate increases as the lattice spacing is reduced. < ipip > at 
the larger lattice spacing L = 4 decays with c = 0.48(5) but at the smaller lattice spacing 
L = 10 decays faster with c = 1.1(1). All fits have x 2 P er degree of freedom w 1 — 2. 

Finally, it should be noted that although the above data is consistent with exponential 
decay, other types of decay can not be excluded. This is mainly due to the fact that since 
the L s = oo result has already being approached, within statistics, for L s = 8 — 10, the 
decay is basically dictated only by the two points L s = 4, 6. More statistics are needed in 
order to clearly establish the type of decay. 



8 Conclusions 

In this paper the properties of Domain Wall Fermions (DWF) were studied in the context 
of the two flavor Lattice Schwinger model. The expectation value of the chiral condensate 
< ipip > was used to probe issues related to restoration of the regularization induced chiral 
symmetry breaking. The expectation value of the relevant t' Hooft vertex < w > was used 
to probe issues related to topology. 
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Dynamical numerical simulations of the full theory were performed. It was found that, 
as expected from perturbative considerations, the restoration of chiral symmetry as a func- 
tion of L s (L s is the size in lattice units of the 2n + 1 direction) at a fixed physical volume 
and lattice spacing, is consistent with exponential decay. In particular, the data is consis- 
tent with a picture where chiral symmetry is restored with a fast exponential decay rate for 
L s up to some value and with a slower exponential decay rate for L s above that value. For 
the range of lattice spacings used in this paper the inflection appeared at L s « 10. Using 
a simple model it was found that the first fast decay can be associated with restoration 
of chiral symmetry in the zero topological sector while the second slower decay can be 
associated with the regions of gauge field configuration space that connect the q = and 
q = ±1 topological sectors. 

The effects of the size of the lattice spacing a to the two decays were studied using both 
analytical arguments and explicit numerical simulations of the full theory. It was found 
that for zero explicit fermion mass the fast decay associated with the zero topological sector 
becomes faster as the lattice spacing is decreased. The vanishing of the chiral condensate is 
consistent with a form e~ cLs with e _c being roughly a linear function of a, but more data at 
smaller lattice spacings are needed before one can be confident that scaling has set in and 
that this is the correct scaling form. The second slower decay also becomes faster as the 
lattice spacing is decreased and it differs less from the faster decay as the lattice spacing 
becomes smaller. For the smallest lattice spacing studied the two exponential decay rates 
differed by « 31 ± 2%. _ 

For small but non zero explicit fermion mass rrif the values of < ipip > and < w > were 
measured. The corresponding L s = oo numbers were calculated by performing numerical 
simulations with the Overlap formalism. It was found that the L s = oo numbers were 
also approached in a way that is consistent with exponential decay with a rate that became 
faster as the lattice spacing decreased. Furthermore, the larger the fermion mass the sooner 
the L s = oo value was approached and for the fermion masses studied in this paper the 
L s = oo result was already achieved to within a few percent for L s = 4 — 8. If a second 
slower decay does set in for 10 ^ L s , it is unimportant and was lost in the statistical 
noise. Finally, an interesting result was obtained from the measurements of < > and 
< w > vs. rrif. It was found that these measurements are in agreement with the analytical 
predictions of |TjJ. In particular the interesting < ipip >~ mj behavior was reproduced. 

Although all the numerical data are consistent with exponential decay, power law decay 
can be excluded only for the fast decay at the smallest lattice spacing studied. For that 
data the decay is sufficiently fast and the error bars are sufficiently small so that a power 
law fit can be safely excluded since it has a y 2 per degree of freedom « 31. More statistics 
and larger L s are needed in order to be able to clearly distinguish between exponential and 
power law decay for the rest of the data points. 

The next step is to carry out a similar investigation for dynamical QCD. Many of the 
characteristics of DWF found here are sufficiently generic so that one would expect that 
they will also be present in QCD. If it turns out that QCD at the presently accessible 
lattice spacings, volumes and quark masses has similar restoration rates as the ones found 
here, then DWF will indeed provide a powerful fermion discretization method. 
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Figure 1. Observables (for definitions see section along the configuration space trajec- 



tory, eq. |66|, labeled by r. The parameters have values: L = 6, fxl = 3.0, m = 0.9, rrif = 0, 
L. = 14. 
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Figure 2. Observables (for definitions see section [3]) along the configuration space trajec- 



tory, eq. |66|, labeled by r. The parameters have values: L = 6, fxl = 3.0, m = 0.9, rrif = 0, 
L s = 4, 6, 8, 10, 12, 14, oo. At L s = oo P 1 = 0. 
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Figure 3. P\ from fig. 2 vs. L s at different values of r = 0, 0.3, 0.4, 0.55 corresponding 
to diamonds, squares, crosses, octagons. 
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Figure 4. Expectation value of the chiral condensate in a configuration space restricted on 
the trajectory of eq. ^ vs. L s for different gauge couplings. The parameters have values: 
L = 6, m = 0.9, m/ = 0, fil = 3.0,2.5,2.0,1.5, and 1.0 corresponding to diamonds, 
squares, crosses, octagons, and stars. 
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Figure 5. Pi along the configuration space trajectory, eq. BB, labeled by r. The parameters 
have values: L = 6, fil = 3.0, m = 0.9, rrif = 0.1, L s = 4, 6, 8, 10, 12, 14, 16, oo. 
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Figure 6. < ipip >l s — < ipip >oo / < ipip >l s vs. L s for rrif = 0.1. < ipip > is the chiral 



condensate in a configuration space restricted on the trajectory of eq. |66|. The parameters 
have values: L = 6, fil = 3.0, and m = 0.9. 
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small fluctuations of size e to the trivial configuration (U = e tr7T with r randomly distributed 
in — e < r < e). The parameters have values: L = 6, m = 0.9, m/ = 0, e = 0.4, 0.3, 0.2, 0.1 
and 0.01 corresponding to octagons, stars, squares, crosses, and diamonds. 
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Figure 8. < w > jvr? vs. mj for L = 6, m = 0.9, fil = 3.0 using the Overlap. The 
dotted lines are the m/ = result ± the corresponding statistical error. 
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Figure 9. < ipip > jm 1 vs. m/ for L = 6, m = 0.9, fil = 3.0 using the Overlap. The 
fits are to < ipip >= AmJ. Both fits have a \ 2 P er degree of freedom of about one. For 
m/ < 0.1 p = 0.996(3), while for m f > 0.1 p = 0.32(2). 
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scale of the mf = 0.05 and nif = 0.0 graphs. The parameters are L — 6, ill — 3.0, m = 0.9 
and L K = 14. 
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Figure 11. < ipip > /m 7 , eq. [4^, vs. L s for mj = 0, fixed physical volume jjl = 3.0, niQ = 
0.9 and for four different lattice spacings fil/L = /xa, with L = 6, 8, 10, 12 corresponding to 
the lines from top to bottom. 



1.0 



0.8 



0.6 



I 



0.4 



0.1 0.2 
1/L 



0.3 



Figure 12. The exponentiated slopes, e c , of the lines in figure 11 vs. 1/L ~ a. The 
diamonds correspond to the faster rates 6 < L s < 10 while the crosses to the slower ones 
12 < L s < 22. 
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Figure 13. [< ipip > /m 7 ] 3 vs. L s for four different mf. The physical volume is fixed 
at /il = 3.0, m = 0.9 and the lattice spacing is set by L = 6. The fits are to a function 
A + Be~ cLa . The dotted line is the L s = oo result ± the error from figure 9. The cross is 
the coefficient A. 
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Figure 14. < w > jm? vs. L s for four different irif. The physical volume is fixed at 
/il = 3.0, m Q = 0.9 and the lattice spacing is set by L = 6. The fits are to a function 
A + Be~ cLs . The dotted line is the L s = oo result ± the error from figure 8. The cross is 
the coefficient A. 
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Figure 15. [< ipip > /m 7 ] 3 vs. L s for four different lattice spacings set by L at fixed 
rrif = 0.2. The physical volume is fixed at ytl = 3.0 and m = 0.9. The fits are to a 
function A + Be~ cLs . The dotted line is the L s = oo result ± the error. The cross is the 
coefficient A. 
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Figure 16. < w > jvr? vs. L s for four different lattice spacings set by L at fixed mf = 0.2. 
The physical volume is fixed at fil = 3.0 and m$ = 0.9. The fits are to a function A+Be~ cLs . 
The dotted line is the L s = oo result ± the error. The cross is the coefficient A. 
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Figure 17. [< ipip > /m 7 ] 3 vs. L s for two different lattice spacings set by L = 4 (squares) 
and L = 10 (diamonds) at m = 0.9. The physical volume and rrifL are fixed at jil = 3.0 
and rrifL — 2.0. The fits are to A + Be~ cLs . The dotted lines are the L s = oo results ± 
the error. The cross is the coefficient A. 
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Figure 18. < w > /m 2 vs. L s . The parameters are as in figure 17. 
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